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Abstract 

• i-H 

ryj The oscillator, inherently, turns the phase noise of its internal com- 

ponents into frequency noise, which results into a multiplication by l// 2 
C/) in the phase-noise power spectral density. This phenomenon is known as 

the Leeson effect. This report extends the Leeson effect to the analysis 
of amplitude noise. This is done by analyzing the slow- varying complex 
envelope, after freezing the carrier. In the case of amplitude noise, the 
classical analysis based on the frequency-domain transfer function is pos- 
sible only after solving and linearizing the complete differential equation 
^ that describes the oscillator. The theory predicts that AM noise gives an 

^\ additional contribution to phase noise. Beside the detailed description of 

CO the traditional oscillator, based on the resonator governed by a second- 

order differential equation (microwave cavity, quartz oscillators etc.), this 
report is a theoretical framework for the analysis of other oscillators, like 
for example the masers, lasers, and opto-electronic oscillators. 
f^*) This manuscript is intended as a standalone report, and also as com- 

plement to the book E. Rubiola, Phase Noise and Frequency Stability in 
Oscillators, Cambridge University Press, 2008. ISBN 978-0-521-88677-2 
(hardback), 978-0-521-15328-7 (paperback). 
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Notation 

Symbol Meaning 

A amplifier voltage gain (thus, the power gain is A 2 ) 

A(s) normalizcd-amplitude noise, see a(t) 

bi coefficients of the power- law approximation of S v (f), 

b(t) <H- B(s) resonator impulse response. Also h v {t) B v (s), etc. 

£ (s) normalizcd-amplitude noise, see e{t) 

f Fourier frequency, Hz 

f c amplifier corner frequency, Hz 

/l Leeson frequency, Hz 

F amplifier noise figure 

h(t) o H(s) impulse response 

h(t) •<->• H(s) phase or amplitude impulse response, takes subscript ip or a 

hi coefficients of the power-law model of S a (f) or S y (f) 

(in case of ambiguity use hi for S y (f) and h\ for S a (f)) 

j imaginary unit, j 2 = — 1 

k Boltzmann constant, k = 1.381 xlO~ 23 J/K 

C{ } Laplace-transform operator 

single-sideband noise spectrum, dBc/Hz. 

Sf(f) - lS v (f), by definition 

Af(s) gain fluctuation, see r)(t) 

Q resonator quality factor 

R, Ro resistance, load resistance (often, R n = 50 O) 

s Laplace complex variable, s = a + ju; 

s derivative operator 

S a (f) one-sided power spectral density (PSD) of the quantity a 

t time 

T, To absolute temperature, reference temperature To = 290 K 

u(t) Heaviside (step) function, u(t) = J 5{t')dt' 

u(t) <-» U(s) voltage (amplifier input) 

v(t) o V(s) voltage. Also a dimcnsionlcss signal 

v(t) complex-envelope associated to v(t) 

V, Vo dc or peak voltage 

V , V Q (t) phasor associated with an ac signal v(t) 

x(t) generic function 

y(t) generic function 

y(t) fractional frequency fluctuation, y(t) = [v(t) — ^ol/^o 

a(t) -H- A(s) normalizcd-amplitude noise, may take subscript u or v 

P(s) transfer function of the feedback path 

7 amplifier compression parameter (0 < 7 < 1) 

6(t) Dirac delta function 

e(t) -f-> £{s) normalizcd-amplitude noise 

rj(t) "R- Af(s) amplifier gain fluctuation 

k small phase or amplitude step 

v frequency (Hz), used for carriers 

a real part of the Laplace variable s = a + jcu 

t resonator relaxation time 
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ip(t) «-» $(s) 


phase noise 


X 


dissonance, \ = oj/oj n — u> n /u> 


V>(<) -o- *(s) 


phase noise (amplifier) 


u; 


imaginary part of the Laplace variable s = o + ju> 


(j 


angular frequency, carrier or Fourier 


UJQ 


oscillator angular frequency 




resonator natural angular frequency 




resonator free-decay angular pseudo-frequency 


f2 


detuning angular frequency, = ujq — w n 


Subscript 


Meaning 





oscillator carrier, as in w , Pq, Vq, etc. 


i 


input, as m Vi(t), <pi\t), ^i(s) 


i 


electrical current, as m the shot noise bi(u>) = 2qi 


n 


resonator natural frequency (u n , u n ) 


o 


output, as in v (t\ tp (t), <& (s) 


P 


resonator free-decay pseudo-frequency (ui p , v p ) 


Symbol 


Meaning 


< > 


mean. Also < >n mean of N values 


— 


time average, as in x 




transform inverse-transform pair, as in x(t) «-> X(s) 


* 


1 i • * / i \ T/i\ T T / \ TT / \ T T f \ 

convolution, as in v(t) = ti(t) * u o V(s) = H(s)U (s) 




asymptotically equal 


O 


zero of a function (complex plane, Bode plot, or spectrum) 


X 


pole of a function (complex plane, Bode plot, or spectrum) 


Acronym 


Meaning 


AM 


Amplitude Modulation (often 'AM noise') 


CAD 


Computer- Aid Design (software) 


FM 


Frequency Modulation (often 'FM noise') 


PM 


Phase Modulation (often 'PM noise') 


PSD 


(single-side) Power Spectral Density 


RF 


Radio Frequency 
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1 Introduction 

The oscillator noise, which in the absence of environmental or aging effect is 
cyclostationary, is best described as a baseband process, after freezing the pe- 
riodic oscillation. The polar-coordinate representation of the limit cycle splits 
the model of the oscillator into two subsystems, in which all signals are the 
amplitude and the phase of the main system, respectively. Putting things sim- 
ply, these two subsystems are (almost) decoupled and all the nonlinearity goes 
to amplitude. This occurs because amplitude nonlinearity is necessary for the 
oscillation to be stationary. Conversely the phase, which ultimately is time, 
cannot be stretched. 

The baseband equivalent of a resonator, either for phase or amplitude is a 
lowpass filter whose time constant is equal to the resonator's relaxation time. 
Hence, the phase model of the oscillator consists of an amplifier of gain exactly 
equal to one and the lowpass filter in the feedback path, as extensively discussed 
in [1]. The amplitude model is a nonlinear amplifier, whose gain is equal to 
one at the stationary amplitude and decreases with power, and the lowpass 
filter in the feedback path. In the baseband representation both AM and PM 
perturbations map into additive noise, even in the case of flicker and other 
parametric noises. This model gives a new perspective on the classical van der 
Pol oscillator. 

The elementary theory of nonlinear differential equations tells us that nonlin- 
earity stretches the feedback time constant. Asymptotically, the time constant 
is split into two constants, one at the oscillator startup and one in stationary 
conditions. If the gain varies linearly with amplitude, which is always true for 
small perturbations, the oscillator can be solved in closed form. 

After the pioneering work of D. B. Leeson [2j, a number of different analyses 
has been published. Sauvage derived a formula that has the same behavior of 
the Leeson formula using the autocorrelation functions [3 . Hajimiri and Lee 
[U 03 |B] proposed a model based on the "impulse-sensitivity function" (ISF), 
which emphasizes that the impulse has the largest effect on phase noise if it 
occurs at the zero-crossing of the carrier. This model, mainly oriented to the 
description of phase noise in CMOS circuits, is extended in [7]. Demir & al. 
proposed a theory based on the stochastic calculus [5] , in which they introduce 
a decomposition of phase and amplitude noise through a projection onto the 
periodic time- varying eigenvectors (the Floquet eigenvectors), by which they 
analyze the oscillator phase noise as a stochastic-diffusion problem. This theory 
was extended to the case of 1// noise [9]. Other articles are mainly oriented 
to the microwave oscillators [TUl El [H] . Demir inspired a work on the phase 
noise in opto-electronic oscillators [13 . Some of the articles cited make use of 
sophisticated mathematics, as compared to our simple methods. All give little 
or no attention to amplitude noise. 

It is well known that the instability of the resonator natural frequency con- 
tributes to the oscillator fluctuations, which in some cases turns out to be the 
most important source of frequency fluctuations. Almost nothing is known 
about the amplitude fluctuation of the resonator. That said, the resonator 
instability is not considered here. This article stands upon our earlier works 
CQ and [Hj. The latter is mainly oriented to the ultra-stable quartz oscillator. 
Here, we present an unified approach to AM and PM noise in oscillators by ana- 
lyzing the mechanism with which the noise of the oscillator internal components 
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is transferred to the output. This theory turns out to be particularly suitable 
to high stability oscillators, based on high quality-factor quartz resonators, mi- 
crowave whispering gallery resonators, etc. This work is easily extended to the 
delay-line oscillator, including the opto-electronic oscillator [15], I16j. 

2 Basics 

2.1 Phase noise 

Phase noise is a well established subject, clearly explained in classical references, 
among which we prefer [TTl [18l [T9l [20] and [21] vol. 1, chap. 2]. The reader may 
also find useful [1, chap. 1]. 

The quasi-perfect sinusoidal signal of angular frequencjQwoj of random am- 
plitude fluctuation a(t), and of random phase fluctuation ip(t) is 

v(t) = [1 + ci(t)] cos [<jj t + (/?(£)] (clock signal) . (1) 

We may need that \ct(t)\ <C 1 and |<,f(i)| -^i 1 or \tp(t)\ <C 1 during the mea- 
surement. The phase noise is generally measured as the average PSD (power 
spectral density) 

SM) = (Mjf)\ 2 ) m (avg, m spectra), (2) 

The uppercase denotes the Fourier transform, so <p(t) o <&(jf) form a trans- 
form inverse-transform pair. In experimental science, the single-sided PSD is 
preferred to the two-sided PSD because the negative frequencies are redundant. 
It has been found that the power-law model describes accurately the oscillator 
phase noise 

o 

SM) = E b ^ n (p° wer law ) ( 3 ) 

n=— 4 

coefficient noise type 

&_4 frequency random walk 

fo_3 flicker of frequency 

fo_2 white frequency noise, or phase random walk 

&_i flicker of phase 

&o white phase noise 

The power law relies on the fact that white (/°) and flicker (1//) noises exist 
per-se, and that phase integration is present in oscillators, which multiplies the 
spectrum xl// 2 . This will be discussed extensively in Section [SJ Additional 
terms with n < —4 show up at very low frequency. 

The two-port components are also described by the power-law ([3]). Yet at 
very low frequency the dominant terms cannot be steeper than / , otherwise 
the input-output delay would diverge. Nonetheless, these steeper terms may 
show up in some regions of the spectrum, due to a variety of phenomena. 

1 ~We use interchangeably ui as a shorthand for 2-kv for the carrier frequency, and as a short- 
hand for 2nf for the offset (Fourier) frequency, making the meaning clear with appropriate 
subscripts when needed but omitting the word 'angular.' 
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2.2 Amplitude noise 

Amplitude noise, far less studied than phase noise, can be described in very 
similar way. The reader may find useful the report [22j . which describes in 
depth the experimental aspects. 

The amplitude noise is generally measured as the average PSD 

S a (f) = (\A(jf)\ 2 ) m (avg, m spectra) , (4) 

where a(t) ■<-> A(juj) form a transform-inverse-transform pair, and described in 
terms of power law 

o 

S*(f) = ]T Kf n (power law) (5) 

n=— 4 

coefficient noise type 

h-3 (integrated flicker) 
h-2 amplitude random walk 
h-i flicker of amplitude 
ho white amplitude noise 

In the current literature the coefficients hi are used to model S y (f), i.e., the 
PSD of the fractional- frequency fluctuation y(t). We use the same coefficients for 
simplicity, because the relationships between spectrum and two-sample (Allan) 
variance are the same. 

At very low frequency, the amplitude noise cannot be steeper than f^ 1 oth- 
erwise the amplitude would diverge. This applies to both two-port components 
and oscillators. Yet, steeper terms can show up in some regions of the spectrum. 

2.3 Oscillator fundamentals 

The simplest form of oscillator is a resonator with an amplifier of gain A in closed 
loop that compensates for the resonator loss[^]l//3. Stationary oscillation takes 
place at the frequency ujq that verify A/3 = 1. This is known as the Barkhausen 
condition. The actual oscillator can be represented with the scheme of Fig. [T] 
which includes a gain compression mechanism and noise. For our purposes the 
noise is represented in polar coordinates as amplitude noise and phase noise. 
The gain compression is necessary for the amplitude not to decay or diverge. 
We assume that A is independent of frequency, at least in a range sufficiently 
larger than the resonator bandwidth. For the sake of simplicity we normalize 
the loop elements so that A = 1 and /3 = 1 at the oscillation frequency u = ujq 
and at the nominal output amplitude v = 1. 

The resonator has narrow bandwidth, hence it eliminates all the harmonics 
multiple of luq generated by the amplifier nonlinearity. Though the harmonics 
can be present at the amplifier output, where the signal can be distorted, they 
do not participate to the regeneration process that entertains the oscillation. 
Hence, the only practical effect of the nonlinearity on the loop dynamics is to 
reduce the gain at the fundamental frequency wo- Therefore, the quasi-sinusoidal 
approximation can be used. 



2 Since the quantity /3 is the resonator gain, 1//3 is the loss. 
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low-pass 



*- b(t) ~ B(s) 
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a u (t) ~ 5^(3) 



E(t) 
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gain fluctuat. 
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— b(t) «• B(s) 



■J saturation 
A 



» L.."j 



Figure 1: Feedback oscillator and its decomposition in PM and AM models. 



Assuming that, as it occurs in practice, the resonator relaxation time r is 
larger than I/wq by a factor of at least 10 2 , the oscillator behavior can be 
mathematically described in terms of the slow-varying complex envelope, as 
amplitude and phase were decoupled from the oscillation. In this representation 
the oscillator splits into two subsystems, one for phase and one for amplitude, 
as shown in Fig. [l] Since phase represents time, which cannot be stretcheoj^] all 
the non-linearity goes in the amplitude subsystem. 

The main advantage of the slow- varying envelope representation is that am- 
plitude noise and phase noise can be represented as additive noise phenomena, 
regardless of the physical origin. This eliminates the difficulty of flicker noise 

3 This is no longer true in extreme nonlinear oscillators, like the femtosecond laser, which 
are out of our scope. 
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™" 1 U 

Figure 2: Most common types of gain saturation. The quantities u and v are 
the rms amplitude at the carrier frequency. 

and other parametric processes. The formalism is simple and tightly connected 
to the experimentally observable quantities. 

We are interested in the mechanism that governs the noise propagation of the 
internal components to the oscillator output. Virtually all oscillators are stable 
enough for the noise to be a small perturbation to the stationary oscillations, 
and consequently for a linear model to be accurate for any practical purpose. 
Linearization gives access to the Laplace-Heaviside formalism. The respons^] 
y(t) to the input x{t) is therefore given by 

y(t) = x(t) * h(t) o Y{s) = X(s)H{s) 

where h(t) is the impulse response, i.e., the response to the Dirac 5(t) function, 
H(s) is the transfer function, the symbol '*' is the convolution operator, the 
double arrow '-fV stands for Laplace transform inverse-transform pair, and s = 
a + juj is the Laplace complex variable. Given the input power spectral density 
S x (f), the output power spectral density is given by 

S y (f) = \H(jf)\ 2 S x (f) . 

The application of this idea to the oscillator rises some difficulties, which will 
be solved in the next Sections. 

3 Amplifier saturation and noise 
3.1 Gain compression 

In large signal conditions, all amplifiers have some kind of nonlinearity that lim- 
its the maximum output power. Neglecting the band limitation, when a sinu- 
soidal signal Urf(t) = U\ cos(wo£) is present at the input of an amplifier, the sat- 
urated output can be written as the Fourier series V T i(t) = Xm=i cos(na;oi). 



4 Here x(t) and y(t ) are generic functions of time, thus not the phase time and the fractional 
frequency fluctuation commonly used in the oscillator literature. 
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Figure 3: Most common types of gain saturation. The quantities u and v are 
the rms amplitude at the carrier frequency. 

The n — 1 term is the fundamental and the n > 1 terms are the harmonics gen- 
erated by nonlinearity. The effect of the band limitation is to change (in most 
cases to reduce) the amplitudes V n and to introduce a phase in each sinusoidal 
term. In a linear amplifier only the fundamental is present at the output, thus 
V n = 0, Vn > 2. 

In the case of the oscillator, the resonator allows only the fundamental to be 
fed back to the input, for the harmonics can be neglected. Hence, the oscillation 
amplitude is described using the slow-varying signals u and v instead of the 
instantaneous peak amplitudes U\ and V\ . Let us define the amplifier gain as 

v 

A = — (definition of A) , 
u 

which of course is equivalent to A = V\/U\. The gain A should not be mistaken 
for the differential gain dv/du. 

Figure [2] shows the gain-saturation types most frequently encountered and 
described underneath. The small-signal gain is denoted with Aq and the gain 
at the oscillator nominal amplitude u = 1 is denoted with A\ . Figure [2] is 
normalized for A\ — \. Around u = 1 the gain can be linearized as 

A = A x [1 - j(u ~ 1)] whith < 7 < 1 , 

which rewrites as 

A = 1 - j(u - 1) whith < 7 < 1 (6) 

after normalizing for A± = 1. 

The slope —7 deserves some comments. The condition 7 > is obvi- 
ous because the gain A(u) must decrease monotonically (Fig. [2]). This is the 
amplitude-stabilization mechanism. A second obvious condition is that in the 
regular-operation region (i.e., around u = 1) the output v(u) must increase 
monotonically. We show that this second condition is equivalent to 7 < 1 by 
substituting ^ in v = Au 

v = - 1U 2 + (1 - 7) (Fig. § , 

The latter is the 'cap' parabola shown in Fig. [3j For v(u) to increase mono- 
tonically at it = 1 it is necessary that the point u = 1 is on the left-side of the 
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maximum, thus 



whose solution is 7 < 1. 

The property that v(u) increases monotonically holds for virtually all am- 
plifiers. Only a few exceptions are found, the most remarkable of which is a 
channel that includes a Mach Zehnder electro-optic modulator. In such cases 
the Barkhausen A(3 = 1 may be met (at least) at two different amplitude levels, 
the first with 7 < 1 and the second with 7 > 1. In one of such cases, it has 
been mathematically proved and experimentally observed that the oscillation 
amplitude flips between these two levels, producing an amplitude oscillation at 
half the frequency determined by the loop roundtrip time [23] . 

3.2 Gain saturation in real amplifiers 

3.2.1 Quadratic (van der Pol) amplifier 

In the classical van der Pol oscillator |24j . the amplifier input-output function 
is defined as y — a\x — 03a; 3 , with a\ > and 03 > 0. In mathematical treatises 
the coefficients a\ and 03 are sometimes set to one. Feeding the signal U r f(t) = 
U± cos(woi) i n such amplifier and taking only the fundamental frequency, the 
output is V r f(t) = Ui [eti — asjUi] cos(cjQt). Accordingly, the gain becomes 
A = ai — a^jUi , which is a 'cap' parabola. 

3.2.2 Hard-clipping amplifier 

In small-signal condition the gain is Aq, independent of the signal level. In- 
creasing the input level the output is clipped when it hits a threshold, where 
the sinusoid progressively turns into a square wave. The asymptotic amplitude 
of the fundamental is 4/tt (2.1 dB) higher than the threshold. This behavior is 
often encountered in amplifiers linearized by a strong feedback, as most circuits 
based on operational amplifiers. Of course the feedback is no longer effective 
when the output is expected to exceed the supply voltage. 

3.2.3 Soft-clipping amplifier 

With moderate feedback, the output clipping starts gradually when the output 
approaches the dynamic-range boundary. This behavior is typical of microwave 
amplifiers. The knee of the gain curve occurs approximately at the 1 dB com- 
pression power. 

3.2.4 Linear-compression amplifier 

The gain law A = 1 — 7(1* — 1) holds in the whole dynamic range. However this 
model may seem a mere academic exercise, it provides useful results in a simple 
and compact form. 

3.3 Amplitude and phase noise 

The contents of this Section is extensively discussed in [25] , and briefly summa- 
rized here. 
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v n (t) = v x (t) cos(a) t) - v y (t) sin(ai t) 



Figure 4: Additive noise. 



3.3.1 Additive noise 

Let us consider a quasi-perfect device that adds a noise term v n (t) to the sinu- 
soidal input signal, as shown in Fig. [4] Assuming that the device gain is equal 
to one, the output signal is 

v(t) = Vq coswoi + v n (t) , 

or equivalently 

v(t) = Vq cos LUot + v x (t) cos (jj^t — v y (t) sin ujQt . (7) 

The random variables v x (t) and v y (t), called in-phase and quadrature compo- 
nent of noise, represent the noise v n (t) in the bandwidth of interest. 

Though the Cartesian representation ^ is the closest to the physics of 
additive noise, polar coordinates can also be used 

v(t) =V [1 + «(«)] cos [wot + ¥>(*)] , (8) 

where in low-noise conditions it holds that 

o(i) = -— and ^(t) - ^— . 

The most relevant feature of the additive noise is that all the statistical 
properties of v n (t), thus of v x (t) and v y (t), are not affected by the input signal. 
There follow some relevant properties 

1. Referred to the input, there is an equal amount of AM and PM noise. 
Yet, AM/PM asymmetry can show up at the output if the amplitude 
non-linearity compresses the AM noise. 

2. AM and PM noise are statistically independent. 

3. The shape of the noise spectrum is independent of the carrier frequency 
u)q. Therefore the noise spectrum cannot have a term 1//, l// 2 etc., 
centered at an arbitrary carrier frequency uiq. 

4. The AM noise and the PM noise scale down with the carrier power. 

The additive noise is generally white, though it can have bumps due to device 
internal structure and it rolls off out of the bandwidth. In the case of thermal 
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Figure 5: Parametric noise. 



noise, including the noise figure F (defined only at the temperature T = 290 
K), the noise PSD of a generic white- noise process n(t) is 

S n (f) = FkT . (9) 
In polar coordinates the noise PSD is S a (f) = ho and S v (f) = bo, with 



hn = 



FkTp 
Po 



and 



bo = 



FkT 
Po ' 



(10) 



where Po is the carrier power. 

In the case of cascaded amplifiers, the Friis formula applies, by which the 
noise contribution of each stage is divided by the gain of the preceding stages 



h 
b 



kTo 
Po 

kTo 
Po 



Fx 



Fx 



Fo 



A\ A\A\ 



Fa 



Fo 



A\ A\A\ 



2 A 2 



(Friis, AM noise) (11) 
(Friis, PM noise) (12) 



where Po is the carrier power. 



3.3.2 Parametric noise 

The parametric noise originates from a near-dc process that modulates the car- 
rier, as shown in Fig. [5] Accordingly, the polar-coordinate representation ^ 
is the closest to the physical mechanism. The most important parametric ran- 
dom phenomenon is flicker noise, whose PSD is proportional to 1/f over sev- 
eral decades. Other types of parametric noise, with PSD proportional to \ j p , 
i = 2, 3, . . ., can only exist in a limited frequency region. For example, l// 5 
noise in the region between 1 mHz and 1 Hz has been observed as the phase noise 
of radio-frequency amplifiers, and also as the offset fluctuation of operational 
amplifiers. If these high-slope phenomena would be allowed to span over many 
decades at low frequencies, the amplitude or the group delay would diverge in 
the long run, which does not fit the experience about two-port devices. 

As a realistic approximation, one can assume that the near-dc process and 
the modulation efficiency are independent of the carrier power P , hence the 
the statistical properties of a(t) and ip(t) tend to be constant in a wide range 
of power 

h-i — C\ and 6_i = C*2 constant, independent of Pq • (13) 
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Nonetheless, a too large input power may affect the dc working point, and in 
turn the amount of parametric noise. 

That the parametric noise is independent of Pq has the amazing consequence 
that the noise of a multistage amplifier is the sum of the individual contributions 

S a (f) = [S a (f)] 1 + [S a (f)] 3 + ... (14) 
SM) = [^(/)] 1 + [5 v (/)] 2 + ... , (15) 

independently of the order of the single stages in the chain. 

Generally, a parametric process affects both amplitude and phase with sep- 
arate coefficients, as the dashed noise generator in Fig. [5] does. This introduces 
some correlation between AM and PM noise. Evidence of this statement is 
provided by the following examples. 

• In a bipolar transistor, a noise source may affect the thickness of the base 
region. Such a process modulates simultaneously the gain (AM noise) and 
the BE BC capacitances (PM noise), which turns in fully-correlated AM 
and PM noise. 

• In a laser medium, the pump power affects the partition between excited 
atoms and ground-state atoms. Two noise phenomena are simultaneously 
driven by the power fluctuation of the pump. The first and more obvious 
phenomenon is the fluctuation of gain and of saturation power, which 
shows up as AM noise — referred to as RIN in the jargon of laser optics. 
The second phenomenon results from the fact that the contribution of an 
atom to the refraction index changes if the atom is excited. This produces 
phase noise inside the loop, thus frequency noise in the laser beam. 

• The third example is provided by the fluctuation of cathodic emission in 
vacuum tubes, like triodes, klystrons, magnetrons, TWTs, etc. Beside the 
obvious effect on gain, the electron emission impacts on the space charge, 
and in turn on the capacitance seen by the signal. 

In all the above examples, a single phenomenon yields fully correlated amplitude 
and phase noise. 



3.3.3 Phase and amplitude noise spectrum 

It has been seen that the AM noise and the PM noise spectra are 

S a (f) = ho + ^~ + ... h = ] -^ ] - h-i=d (constant) (16) 
J -H) 

S v {f)=bo+ h —^ + ... b Q = ^- 6_! = C 2 (constant) (17) 

J M) 

An example of phase noise spectrum is shown in Figure [6] This Figure em- 
phasizes the fact that the flicker noise is constant and that the white (additive) 
noise scales down as the power increases. The obvious consequence is that the 
corner frequency f c also scales with power. A common mistake found in CAD 
software is that the flicker is described by a fixed corner frequency, independent 
of power. The reader is strongly encouraged to check before trusting a CAD 
program. 
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f Q = / FkT Q )P Q depends on P n 



b = FKT Q /P Q 
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Figure 6: Amplifier phase (or amplitude) noise power-spectrum density. 



A ' 1 1 oscillator operation I 



fluctuating 
A = 1-y(u-1)+r| 




1 u 

Figure 7: Parametric fluctuation of the amplifier gain. 



3.3.4 Gain fluctuation 

Modeling the oscillator in the frequency domain, the gain is A — A\[l — 'y(u—l)] 
around the oscillation point u = 1. Introducing a slow fluctuation, A turns into 
the slow varying function of time 

A(t) = A x [1 - 7 (« - 1)] [1 + r)(t)] 
approximated as 

A(t)^A 1 [l- 1 {u-l)+n{t))e^ (Fig.[7f, (18) 

where 

r](t) O A/"(s) (amplitude fluctuation, i.e., AM noise) 
ip(t) O \l/(s) (phase fluctuation, i.e., PM noise) 

are the amplitude an phase gain fluctuations, respectively. Figure [7] shows the 
combined effect of the gain amplitude fluctuation and compression. 

Flicker noise, which results from a parametric effect, impacts directly on the 
gain. It can be described by 



[Sri(f)] flicker 



h-i 



and 



[-%/>(/) J flicker 



h : 



(constant vs. P ) 
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[1+£(t)] COS[U) t+ip(t)] 


resonator 


[1+«(t)] cos[w t+cp(t)] 


► 




— \ — r 



8(t) 



[1+£(t)] COS[U) t+ip(t)] 



resonator 

Q, t« n 



[1+<x(t)]cos[w t+cp(t)] 

t — t 



set a small phase or amplitude step k at t=0, and linearize for k->0 



cos(u> n t) t = resonator 




COS(U)gt+K) 



[1+(x(t)]cos[u) t+cp(t)] 

\ 



k Jb w (t) dt 



k Jb 0[(p (t) dt 



Figure 8: AM and PM response of a resonator. 



Additive noise, albeit of quite different origin, can still be seen as a gain 
fluctuation because it affects the input /output relationship. Hence 



ISM)} 



additive 



additive 



FkT 



(constant vs. f) 



4 The resonator 

The resonator in actual load condition^ is governed, or locally well approxi- 
mated by the differential equation 



+ -^v + u 2 n v = L{vi(t)} 



(19) 



5 Whoever has worked seriously in the field of oscillators, may have in mind three sets 
of parameters like l uj n and l Q.' These sets refer (1) to the unloaded resonator, which is 
a mathematical abstraction not accessible to the physical experiment; (2) to the resonator 
loaded by the measurement test set, from which the unloaded parameters are estimated; and 
(3) to the resonator loaded by the oscillator circuit. The external circuit, either the test set or 
the resonator, increases the dissipation and affects the natural frequency. That said, it is to 
be made clear that here the resonator is always loaded by the oscillator circuit, and therefore 
that there is no point in discussing the other conditions. On the other hand, the process of 
getting £j n and Q from experimental data may be tricky or difficult. We skip this discussion 
because it depends on the specific resonator and oscillator, while we aim at a general theory. 
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where u) n is the natural frequency, Q is the quality factor, i>,(t) is the external 
force, and L is an operator. The most interesting form of the force term is 
L{«i(£)} = ^Q-ii(t) because it is homogeneous with the dissipative term. This 
occurs with the series (parallel) RLC resonator driven by a voltage (current) 
source, and in other relevant cases. Accordingly, (|19[) becomes 



v + + uo 2 n v = -^Vi(t) . (20) 

Using the Laplace transform, the resonator transfer function /3(s) = V /Vi is 

«•>- < 21 » 



Equations ( 20 ) and ( 21 ) are normalized for the resonator to respond to a sinusoid 
at the exact resonant frequency ui n with a sinusoid of the same frequency, phase 
and amplitude. 

We analyze the impulse response of the resonator phase and amplitude in 
stationary-oscillation conditions. The phase response is the response to a pertur- 
bation S(t) in the argument of the driving signal, as shown in Fig. [8] Similarly, 
the amplitude response is the response to a perturbation S(t) in the amplitude 
of the driving signal. In general literature the impulse response is denoted with 
h(t), and its Laplace transform with H(s). Since we use h(t) «-» H(s) for the 
oscillator response, the phase or amplitude impulse response of the resonator is 
denoted with b(t) <H- B(s). It turns out that the resonator response is the same 
for amplitude and phase. 

In our analysis we replace the impulse 5(t) with a small phase or amplitude 
step Ku(t), where u(t) is the Heaviside function 

u(f) = / 5{t) dt = 




and we linearize for k — > 0. Then we use the general property of linear systems 
that the response to u(t) is J b(t) dt. Notice that u(t) can be seen as a switch 
that changes state from off to on at t — 0; and that u(— t) switches from on to 
off at t = 0. 



4.1 Sinusoidal transients 
4.1.1 Switch-off transient 

Let us consider the resonator driven by the signal 

Vi(t) = — cos(uj Q t - 9) (probe) (22) 

Po 

where (3q and 9 are chosen for the asymptotic output to be v (t) — cos(ojo^) for 
t — > oo, i.e., amplitude is equal one and phase equal zero in the general case 
loq 7^ oj n . If the probe signal Vi(t) is switched off at the time t — 0, the output 
is 

v (t) = cos(uj p t) e" 1 ^ t>0, (response) (23) 
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where 



2Q 

t = — (relaxation time) (24) 



is the resonator relaxation time and 



uj p = u) n \/l- -^2 (free-decay pseudo-frequency) (25) 



is the free-decay pseudo-frequency. For Q > 1, we can approximate lo p ~ 
ui n . Thisis justified by the fact that the phase error £ accumulated during the 
relaxation time r is 

C = (Wn - u p )t = — . 

This is seen by replacing r = ^ and cj p = w n y/\ — 1/4Q 2 in £, and by expand- 
ing in series truncated at the first order. 

4.1.2 Switch-on transient 

The response to a sinusoid switched on at the time t = takes the general form 

v (t) = £/ cos(ojpi) + 3§sin(ujpt) e~~ + 'rf cos(ojot) + ^sin(wot) t > , 

where ^, and £F are constants determined by the boundary conditions. 

We use the probe signal (22). This yields immediately 'tf = 1 and 9 = 0. 
The constants ^ and ^ are found by assessing the continuity of v (t) at t = 0, 
which gives ^ = — 1 and SS — 0. Approximating u; p ~ w„ for Q 1, the 
output is 

y (t) = - cos(u;„i) e~ t/r + cos(w i) i > G . (26) 
Similarly, using a probe signal Vi(t) — -i^ sm(ujot — 6), the output is 

v (t) = - sm(u n t) e~ t/T + sin(w t) t > (27) 

4.2 Impulse response at the exact natural frequency 

When the resonator is used at the exact natural frequency, it holds that ujq = u n , 
fa = 1, and = 0. 

A phase step k at £ = is described as the probe signal 

Vi(t) = cos(woi) u(— t) + cos(w £ + «) u(i) , 

v v ' " v ' 

switched off at t — switched on at £ = 



By virtue of linearity, the response is the sum of (23) plus (26), that is, 

v a (t) = cos(wpi) e~* /r + cos(cj p < + «)[l - e~ t/T ] t>0. (28) 

Expanding and using the approximations cos(k) ~ 1 and sin(re) ~ k for n — > 0, 
and ujp ~ uj n for large Q, thus ui p ~ cjo, we get 

^o(i) = cos(wq^) — Ksin(cjoO[ i — e~*/ r ] i > , 
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Resonator response in the phase space 



time, t/tau 

I i i i f i i i 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Figure 9: Resonator response to the step and to the impulse. 



This can be seen as a slowly varying phasor V Q (i) = {l + jft [l 
whose angle 



=.-*/ 



1>. 



-t/r] 



t > 



is the response to nu(t). After deleting k and differentiating, we obtain the 
impulse response b(t) = 



-t/r 



An amplitude step k at t — is described as the probe signal 
Vi(t) = cos(wot) u(— t) + (1 + k) cos(woi) u(t) , 



switched off at t — 



switched on at t — 



Once again the response is the sum of (|23|) plus ( 26 1 



v {t) = cos(wpt) e~ t/r + (1 + k) cos(w p i + k) [l - e"* /r ] t>0. (29) 

Expanding under the same approximations as above, i.e., cos(k) ~ 1 and 
sin(/t) ~ k for k — > 0, and w p ~ w„ for large Q, and w p ~ U)q, we get 

u (t) = cos(w <) + k[1 - e"'/ T ] cos(u! Q t) t > . 

This is a slowly varying phasor V (i) = ^ |1 + k[1 — e _t / T ] } , whose ampli- 
tude swing 



V (*)- V„(0)~k[1 



-t/r] 



f > 



is the response to nu(t). After deleting k and differentiating, we obtain the 
impulse response b(t) = ^e~*/ T , the same already found for the phase impulse. 

In conclusion, the response to either a phase impulse or to an amplitude 
impulse is (Fig. [9]) 



b(t) = -e- 
r 



<-> B(s) = 



1/r 



1/r 



(30) 
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£(t) ~ Sis) 



i|j(t) <-> Y(s) 
•- 



file: sle-AM-PM-cpl-rt 





B acp 












B cpcx 









cpcp 



a v (t) <-> ^(s) 



Cp(t) ^ 0(s) 

• 



Figure 10: Detuning the resonator results in coupling AM to PM. 



Eq. (30) is that of a simple RC low-pass filter, which we will use in all block 
diagrams. The inverse of r is known as the Leeson (cutoff) frequency of the 
resonator 



1 u n 



h = 



1 _ Vn_ 

2-kt ~ 2Q 



Finally, it is to be remarked that ( 28 1 contains no amplitude terms at first 



order, and that (29) contains no phase terms at first order. This means that at 



the exact resonant frequency there is no phase-amplitude coupling. 



4.3 Off-resonance impulse response 

In this Section we analyze the impulse response of the resonator when the carrier 
frequency is luq ^ cu n , with an offset 



W _ W >7 



(detuning) . 



An amplitude perturbation e(t) in the resonator driving signal results in an 
amplitude fluctuation ot{t) — b act (t) * s(t) plus a phase fluctuation ip (t) — 
b va (t) * Similarly, the resonator responds to a phase perturbation tp(t) 

with a phase fluctuation tp(t) — b V(p (t) * tp{t) phis an amplitude fluctuation 
a(t) = b av (t) * tp(t). This is written in matrix form as 



baa 
b[^>a 



Jiptp 



B, 
B 



ff. 



and shown in Fig. 



response is (Fig. 11) 



10| In the following sections we will prove that the step 



[b](t) dt = 



sin(fft) e~*/ r 
[1 - cos(fii)] e"*/ T 



[1 - cos(m)] e~'/ T 
sin(fit) e-*/ r 



(31) 



and that the impulse response is (Fig. 12) 
[b] (t) 



( Q sin Qt 



-t/r 



(— Q cos Oi + i sin Qt) . 



-t/r 



[-Q cos tit + i sin Oi) e~ t/r 
(fi sin D,t + \ cos ttt) e~*/ T 



(32) 
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-0.5 



Step response of the detuned resonator 




time , t/ljau t 
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Figure 11: Step response of the resonator off the natural frequency. 



Impulse response of the detuned resonator 



Parameter: detuning frequency 




time , t/tjau t 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Figure 12: Impulse response of the resonator off the natural frequency. 



[B] (s) = 



1 



s + 



T S z 



-a + - 
-fis 



Q 2 



1 



s 2 + h - 



■ + Q 2 t s 2 + f s + 4r - 

The resonator response has diagonal symmetry 

b aot (t) =b w (t) «-» B QQ (s) = B w (s) 
b tt¥ >(*) = bp a (t) o B Qy (s) = B ¥Q (s) . 

The proof is given in Sections [Olj [4331 and |4~3~4l 



k + n 2 



(33) 



(34) 
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4.3.1 Response to the phase impulse 

A phase step k at t = is described as the probe signal 

Vi(t) = — cos(uj t — 9) u(—t) + — cos(uj t — 9 + k) u(t) 
Po Po 

S v ' S v ' 

switched off at t — switched on at t — 

= — cos(oj t — 9) u(—t) + i [cos(woi — 0) cos k — sin(w t — 0) sin k] u(t) 

Po Po 

■ cos(oj ^ — 8) u(— i) + -r- [cos(woi — 9) — k sin(u;oi — 0)1 u(i) k«1. 

Po 



A) 



Using (23), (26 1 and (27) under the large-Q approximation (uj p — uj n ), the above 
yields the output 

v (t) — cos(w„t) r_ cos(a;„t) + cos(wof)] 

+ K[sin(u;„i) _ sin(u;ot)] (t > 0) 

which simplifies to 

v Q (t) = cos(wo^) ~~ Ksin(a;o<) + nsm(ui n t) e~ l ^ T t > . 

Introducing the detuning frequency f2 = wo~w„, we get sin(w n t) = sin(woi— fM), 
thus sin(w n t) = sin(wot) cos(fM) — cos(wni) sin(fit). Hence, the output signal can 
be rewritten as 

v (t) = cos(cjo^) — Ksin(woi) + Ksin(u)ot) cos(ftt) e~ v ~ k cos(woi) sin(fii) e~ v , 
which simplifies to 

v (t) = cos(w i) 1 — «sm(Q£)e~* / ' T — Ksm(uj Q t) 1 — cos(fi^)e~'/ T . (35) 
Freezing the oscillation uot, the above turns into the slow-varying phasor 



V Q (t) = {l - Ksm(nt)e- t/T + j K [l - cos(m)e~ t/r ] | 



k < 1 



of angle 



arctan = k 



1 - cos(m)e-* /r 



5R{V (<)} 
and amplitude 

|V (i)| = |V o (0)|- K sin(m)e-*/ T 
After deleting k and differentiating, we obtain the impulse response 



i>ipip(t) = f2sin(f2t) H — cos(Oi) 



b a¥ ,(*) = — Qcos(fM) H — sin(fit) 
L r 



phase 
amplitude 



(36) 
(37) 
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4.3.2 Response to the amplitude impulse 

An amplitude step k at t — is described as the probe signal 

Vi{t) — -j- cos(wo* — Q) u (—t) + (1 + K )~B~ cos(o;oi — 9) u(t) 
Po Po 
v v ' v v ' 

switched off at t — switched on at t — 

Using ( 23 ) , ( 26 ) and ( 27 ) , under the approximation uj p = uj n the above yields 
the output 

v a (t) = cos(uj n t) e~ l l T + (1 + k) [- cos(uj n t) e~ t / T + cos(w i)] t > 
= cos(woi) + k cos(ojo^) + k cos(u! n t) e~ l / T . 
Using cos(cj„t) = cos(o;oO cos(f2i) + sin(cjot) sin(ftt), the output is 



v a (t) = cos(w i)|l + 



1 - cos(ftt) e~ t/r \-k sm(uj t) sm{nt)e~ t/T (38) 
Freezing the oscillation u: t, the above turns into the slow-varying phasor 

V Q (t) = {l + «[1 - cos(m)e~ t/T ] + j/«sin(Ot)e-' /T } 

of angle 

S{V (t)} _ .„ : „^„_ t/r 



arctan ■ 



«sin(Qi)e 



5R{V (<)} 

and amplitude swing 

|V (t) - V o (0)| = «[l - cos(fit)e-*/ r 



k < 1 



k < 1 



After deleting k and differentiating, we obtain the impulse response 



baa(i) = risin(fit) H — cos(fit) 
L r 

b VQ (t) = — ficos(Qi) H — sin(fii) 



amplitude 
phase 



(39) 
(40) 



4.3.3 Remark 



Interestingly, the phase noise bandwidth of the resonator increases when the 
resonator is detuned. This is related to the following facts. 



1. When the resonator is detuned, it holds that 



d arg[/3(jo;)] 
duj 



< 



d arg[/3(ju;)] 



(41) 



With lower slope, the oscillator phase noise is higher. 



2. Detuning the resonator, the symmetry of arg[/3(jw)] around the oscillation 
frequency is broken. This explains the frequency overshoot seen in Fig. |12| 
for Q ^ 0. 

3. The step response decays faster when the resonator is detuned. 
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4.3.4 Appendix: Formal derivation of [B(s)] from [b(i)] 

For the sake of completeness, we derive the full expression of [B(s)] from [b(i)], 



that is ( 33 1 from ( 32 1 . Thanks to the symmetry properties ( 34 ) , we only need 
to derive B aQ (s) and B aip (s). 



Using the well known properties 



1 

s + 1/t 
e at f(t) f+ F(s - a) 



we notice that it holds 



,±jm e -t/r „ 



1 



(42) 



B Qa (s) = £{b aa (t)} is found using (42 1, after expanding (39) with the Eulcr 
formulae 

cos(fii) = X - {e m + e~ 3nt ) sin(fit) = ~ (e J ' ot - e~ ]nt ) . (43) 



Thus, 



B aa ( s ) = c{[n^( 



[( jat_ e - j nty 



1 1 



)]e"*/r} 



J2 



s - jfi + i s + in + i 

fi 2 + s/r + l/r 2 

i s + ± + n 2 T 



i 

27 



and finally 



B aa (s) 



T S 2 + + 4r + SI 2 



qed 



(44) 



Similarly, B Q¥> (s) = £ {b ay (t)} is found using (42), after expanding (37) 
with the the Euler formulae (|43|. Thus, 



-fls - uj/t + oj/t 2 

; s + i-jn) ( s + i+jn) ' 



i 

j27 



and finally 



B QV (s) = 



-Os 



•O 2 



qed 



(45) 
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Figure 13: Phase-noise model of the feedback oscillator. 



-e — 35— ► 



-1/r 




Figure 14: Phase-noise transfer function. 



5 The Leeson effect 



Figure 13 shows the phase-noise model of the oscillator. In this figure, all signals 
are the phase fluctuation of the oscillator sinusoidal signal. Here, the resonator 
turns into a lowpass filter of time constant r, as explained in Section]!] A noise- 
free amplifier has a gain exactly equal to one because the amplifier repeats the 
phase of the input signal. The real amplifier introduces the random phase ip(t), 
which in this representation is additive noise, regardless of the physical origin. 
For the sake of simplicity, we put in tp(t) all the phase-noise sources. 
We define the phase-noise transfer function as 



hoo 



*(«) 



Applying the elementary feedback theory to the circuit of Fig. 13 we find 

H( S ) ' 



1 + B(s) 



where B(s) is the resonator transfer function (30), and therefore 

s + 1/r 



H!» = 



(Fig. 141 



(46) 



This is the Leeson effect, by which the oscillator integrates the slow phase fluc- 
tuation, turning it into frequency fluctuation. The phase-noise transfer function 
is plotted in Fig. [14] 
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U — 1 + Ol u 
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low-pass 

— MV^- 

- relaxation time 

t = 2Q/lj 



Figure 15: Amplitude-noise model of the feedback oscillator. 

6 Low-pass model of the oscillator amplitude 

Figure [15] shows the low-pass model that describes the oscillator amplitude. 
Since the gain A depends on amplitude, the Laplace/Heaviside formalism can- 
not be used directly. We first need to linearize the system in the appropriate 
conditions. 



6.1 Differential equation 

Cutting the feedback loop at the amplifier input, we get 

u = e + v 2 , 

where t>2 results from the lowpass filter 

1 



r 2 = - / (yi — v 2 ) dt 
Combining the above equations and replacing v\ = Au and v 2 = u — e, we get 
u- - j (A - l)udt = e + ^ J edt (general IE) , (47) 

thus 

u (A - 1) u = i + -e (general DE) . (48) 



Notice that (|47j)-(48) are general because A is still unspecified. Substituting 
A = 1 - j(u - 1), as in Fig. [2] (|48} becomes 



7, „, 1 
u H — [u — l)u = e H — e 

T T 



with A = 1 - 7(u - 1) 



The system free running is governed by the homogeneous equation 



— (u — lju = 



The solution is 

«(*) 



1 



(solution of ( 50 ) ) 



1 + C*e-T*/ r 

where C is a constant determined by the initial conditions. 



(49) 
(50) 
(51) 
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Figure 16: Oscillator startup. 



6.2 Simplified oscillator model 

A simplified model for the oscillator is obtained by assuming that the linear 
approximation A = 1 — ^{u — 1) holds in the whole amplitude range. One can 
object that this case is only of academic interest because in real amplifiers the 
parameter 7 is constant only in a narrow region around u = 1, as shown in 
Fig. [2] Nonetheless, the general description that follows can be easily adapted 
to practical cases. 

Assuming that A = 1 — 7(1/ — 1), the oscillator amplitude is fully described 
by (51 1, where C is related to the initial u(0) by 

u ^ = TTc - C = W)- 1 - 

Thus, 



«(*) = 7 : \ " • (52) 

u(0) J 

In the absence of a switch-on transient, oscillation starts from noise. Thus, it(0) 
is a small positive quantity, hence l/u(0) — 1 ~ l/u(0) and 

u(t) ~ 1 < u(0) « 1 . (53) 

1 + u(Q) C 

Accordingly the following asymptotic expression hold 

u(t) = u(0) e 7t/T t -> (54) 

u(t) = 1 t -> 00 (55) 



Figure 16 shows the complete oscillation start (52), which saturates to u = 1, 
and the small-signal approximation (54). Figure 17 shows a Spice simulation. 
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SPICE simulation 
Y = 0.1, t = 16/js 
u(0) = 1E-4 ... 1E-2 



3.0ms 



4.0ms 



Figure 17: Simulation of the startup of the model shown in Fig. |15[ assuming 
that the amplifier law A = 1 — 7(1* — 1) holds in the full range. 



6.3 Oscillation soft-start 

Actual oscillators differ from the above simplified model in that the small-signal 
gain follows the law A = 1 — 7(1* — 1) only in the vicinity of u = 1, as shown 
in Fig. [2] In the absence of a general model, we denote the small-signal gain 
with Aq, which is a circuit-specific parameter that we assume to be constant for 
u — )■ 0. Replacing A = A (constant), the homogeneous equation (50) becomes 



The solution is 



u (A - 1) u = 

r 



u{t) = u(0)e {A °- 1)t/T 



(56) 



where it(0) is the small initial condition set by noise. This solution is similar to 



(54), but for the different value of the time constant. 

A number of computer simulations were done independently by R. B. well 
before the approach presented here was developed [27] HE]- This led to the 
preliminary work published in [25]. Figure 18 shows the simulated startup. 
The left-hand side of the envelope, until t w 100 /xs, fits well the theoretical 



prediction (56) 



6.4 Restoring mechanism 

It is interesting to study the amplitude free running when the initial condition 
is set close to the steady state, thus at u(0) = 1 + k with < k <C 1. In 
this conditions the approximation A =1— 7(1*— 1) holds, and therefore the 



amplitude is given by (52 ) with u(0) = 1 + k 



«(*) = ^ ~ HO) = 1 + «) , (57) 

1 e -T f / T 

1 + K 



For small k, this is linearized as 



u(t) = 1 + ne-^/ T , 



(58) 
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Figure 18: Simulated oscillator startup 
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or equivalently 



a u (t) 



ne 



-yt/r 



because a u = u — 1 (Fig. 15). The time constant 

T 

T r = - 

1 



(restoring time) 



(59) 



(60) 



is the oscillator restoring time for amplitude perturbations. Since in virtually 
all amplifiers it holds that < 7 < 1, as widely discussed in Section [XT] it holds 
that r r > r. 



6.5 Amplitude impulse response 

Now, we study the oscillator response to the amplitude impulse e{t) = 5(t) 
occurring when the oscillator is in the steady state it = 1, v = 1. The impulse 
response is the derivative of the step response, linearized for small perturbation. 
Thus, referring to Fig. |15[ we apply at the input the small step 



e(t) = Ku(t) 



< K < 1 



We know from Section |6.4| that the small-signal response is a decaying expo- 
nential e -7 '/ r . Hence the response is completely determined by the initial and 
final values 



u(t) = u(po) + [u(0+) - u(oo)] e~ 7t/T 



t > 



(61) 



Since the perturbation takes time to propagate through the lowpass filter, it 
holds that u(0+) = 1 + k. The final value is obtained by inspection on Fig. 15 
after bypassing the lowpass filter and setting e(t) — k. Thus 



U = t'2 + k 

v 2 = [1 - 7(u - 1)] u 



(adder) 

(amplifier and filter) 



hence 



u = [1 — 7(1* — 1)] u + K 



(t — !> 00) 
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The algebraic solutions 

7 — \/l 2 + 7 + \/l 2 + 

Ul = 2^ U2 = 2^ ' 

and for small re 

K K 

Mi x 112XH — re — > . 

7 7 
It is immediately seen that u\ < and w2 > 0. Hence U2 is the physical solution 



while Mi is discarded. Setting re = 1 (unit step) and using (61), we find the step 
response 



h, ( f } (t) <ft = 1 + i u(i) - f i - 1^ e^ t/r u 
7 \7 / 



(i) . (62) 



Notice that the term '1+' is the steady state before the step is applied. The 
subscript u and the superscript (e) , which refer to the input e and to the output 
u, are introduced to emphasize the difference versus other response functions. 



The impulse response is found by differentiating ( 62 ) 



hV(t) = m + l(j-l) e-^u(t). (63) 

The Laplace transform Hu\s) — £{hu\t)} is found immediately using the 
properties £{S(t)} = 1 and £{e- at } = ^ 



u w r V7 J s+j/t ' 



which simplifies to 



H^(,) = ^f. (64) 

So, the transfer function is completely determined by the roots 

(pole) s p = -7/r f. p = jf L < f L 

(zero) Sz = -1/t f z = f L . 

7 Extension of the Leeson effect to AM noise 

In this Section we study the effect of the parametric fluctuation of the gain by 
introducing the random variable rj(t), as anticipated in Section 3.3.4 and Fig. [7] 



A = 1 - -y(u - 1) + r\ (65) 
r)(t) o Af(s) . 

We linearize the system for low noise, and we search for the transfer functions 
H«>(») = £M ,nd H<->(,) = ^, 

where 

a u (t) A u (s) and a v (t) O A v (s) , 
are the amplitude fluctuations at the amplifier input and output, respectively. 
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Figure 19: Amplitude-noise transfer function (amplifier input). 
7.1 Noise at the amplifier input 



By replacing A = 1 — j(u — 1) + r] in the homogeneous equation (50), we get 



M H (U — ljM = — It . 

T T 

Since u = 1 + <x u , it holds that u = a u and u — 1 = a„, thus 

. . 7 »? 
a u H — a u u = — u . 

T T 

For small fluctuations a u and ?y, we linearize the above using u ~ 1 

• 1 

a u + -a u = -i] . 

T T 

The linearized system can now be described in using the Laplace transforms 

(s+^jA u (s) = ^M(s) , (66) 



which gives the transfer function (Fig. 19 ) 



1/r 



s + 7/V 



(67) 



7.2 Noise at the amplifier output 

We first need to relate a v to a„. This is done by replacing A 
in v = Au 

v = {-j{u - 1) + 1 + rj\ u , 
and by expanding using v — 1 + a v and u — 1 + a u 

1 + a v = 1 + i] - ja u + a u - a u rj - -fa 2 u . 
Neglecting the second-order noise terms a u r\ and a 2 

ctv = (1 - l)a u + V , 



-7(tt— 1) + 1 + »J 
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Figure 20: Amplitude-noise transfer function (oscillator output). 
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SPICE simulation 
Y = 0.1, t = 16/js 
n.=1 E-3 step at t=4 ms 



3.8ms 



4.0ms 



4.2ms 



4.4ms 
Time 



4.6ms 
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Figure 21: Simulation of the step response of the AM model shown in Fig. [15} 
we get 



a v - r\ 



A v (s)-Af{s) 



1 — 7 " 1 — 7 

Then, by replacing the above A u (s) in Eq. (|66| we get 



and finally 



s + 1/t 

s + 7/V 



(68) 



The transfer function H^^s) is shown in Fig. 20 Notice that the case 7 > 1 
(dashed green curve) is not allowed by the condition < 7 < 1 for the amplifier 
gain. 



Figure 21 show a SPICE simulation of the oscillator response to a gain step 
rj = 10~ 3 , assuming that the gain-compression parameter is 7 = 0.1. The rising 
exponential reaches the final value 1 + 77/7 with a time constant T/7, which 



confirms Eq. ( 68 1 



7.3 Predicted spectra 

We calculate the oscillator AM and PM spectra due to the Leeson effect alone. 
The fluctuation of the resonator natural frequency, not accounted in this Section, 
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Figure 22: Phase- noise spectrum in log-log scale. The cutoff frequency generated 
by the zero of the transfer function at s = — 1/t is shown on the horizontal axis. 
The pole at s = cannot be shown on a logarithmic scale. 



may be added afterwards. With the remarkable exception of the laser, virtually 
all practical oscillators are followed by a buffer, which contributes with its own 
noise. Referring to first plot of Fig. 22 (amplifier PM noise), we notice that 
the output buffer has higher flicker and lower white noise than the sustaining 
amplifier. The buffer flicker is higher because the buffer has higher number of 
stages, each of which adds its 1/f phase noise independent of the carrier power 
[Eq. ( 13 )] . Conversely, the buffer white noise is lower because this type of noise 
is additive and the input power is higher at the buffer input [Equations |TT] ) 
( 12 )]. The same is seen on the first plot of Fig. 23 which refers to the amplifier 
AM noise. 



7.3.1 Phase noise 



With reference to Fig. 22 the analysis starts from the sustaining-amplifier noise, 



which shows the flicker corner at / = / c . This noise is turned into the oscillator 



noise by the transfer function H v (f) [Eq. (46)], which is completely described 
by a pole at / = and a zero at Jl = on the Bode plot. 

With a \ow-Q resonator we get the spectrum of the Type 1, where /l > fc- 
At the oscillator output, before buffering, only the slopes 1/f 3 , l// 2 and /° are 
present. The buffer noise is generally not visible because it rises with lower slope 
(1//) on the right hand of the plot. Only in a few special cases, when noise 
special techniques are used to reduce the phase noise of the sustaining amplifier 
[501 |3"T] , thus the gap between the flicker of the buffer and of the sustaining 
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Figure 23: Amplitude-noise spectrum in log-log scale. The cutoff frequencies 
generated by the roots of the transfer function at s = — j/t and s = — 1/r are 
shown on the horizontal axis. 



amplifier is large, some 1 // noise shows up at the buffered output in the region 
around /j,. 

If the resonator Q is higher we get the spectrum of the Type 2, where /j, < f c . 
Before buffering, only the slopes l// 3 , 1// and f° are present. The buffer 1// 
noise shows up because it is higher than the sustaining-amplifier noise and has 
the same slope. 



7.3.2 Amplitude noise 



The amplitude noise (Figure 23 1 is more complex than the phase noise because 
the transfer function H„(/) [Eq. (68 1] is described by two roots, a real pole at 
f = 7T— and a real zero at fr = t^—. Increasing the resonator Q, these roots 

Z7TT J Z7TT ° b ' 

may occur both on the right-hand of f c in the Type-1 spectrum (low Q), one 
on the left hand and the other on the right hand of f c in the Type-2 spectrum 
(medium Q), and both on the left-hand of f c in the Type-3 spectrum (high Q). 
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u = l+u v 

a u (t) <-> a u (s) 




Figure 24: Combined effect of the AM-PM coupling in amplifier and the loop, 
which forces the Barkhausen condition. 



Generally, the buffer l/f noise shows up only in the Type-3 spectrum, in 
the l/f region between f^ and f c . It may also show up in the Type-2 spectrum 
around fi if the gap between the flicker of the buffer and of the sustaining 
amplifier is made large by the use of a noise-degeneration sustaining amplifier. 



8 AM-PM coupling in the amplifier 

We turn our attention to the AM-PM noise coupling mechanism shown in 
Fig. [24] Noise modulates the gain. Yet the Barkhausen condition forces the 
loop gain to be equal to one through the gain-compression mechanism. The 
consequence is that the gain fluctuation is transformed into a fluctuation of the 
oscillation amplitude, and in turn into a fluctuation of the amplifier phase. The 
conclusion is that the phase ip seen by the Leeson effect is the sum of two con- 
tributions, the first comes straight from the amplifier, and the second results 
from the effect on the fluctuating amplitude. The detailed model that follows 



is shown in Figures 24 and 25 and discussed underneath. 

For the sake of simplicity, we assume that the oscillator is tuned at the 
exact natural frequency of the resonator, and we assume that the amplifier is 
perturbed by one dominant source of noise. These hypotheses give a realistic 
picture of the oscillator. 

Denoting with x(t) ■<-> X(s) the near-dc noise process, and introducing the 
modulation efficiency and w^, the amplifier gain is perturbed by a factor 
(1 + w J1 x)e' Wi ' x . Accounting for compression and neglecting the second-order 
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Figure 25: Combined effect of the AM-PM coupling in amplifier. The 
Barkliausen condition turns the gain fluctuation 77 into amplitude fluctuation, 
and in turn to the phase fluctuation -0". The latter adds to the phase fluctuation 
<//. 



terms, the complete expression of the gain is 

A = [1 - 7(1* - 1) + w n x] e? w * x . 

The stationary oscillation is ruled by the Barkhausen condition Af3 = 1. With 
the normalization /3(wo) = 1, this implies that |^4| = 1 . There follows that the 
instantaneous gain fluctuation rj(t) cannot increase Instead, n(t) causes 
the oscillation amplitude u(t) to change from I to I + r/(t)/j, as shown in the 
upper plot of Fig. |25| In this condition the amplifier introduces a phase term 
ip" oc 77/7, which adds to the 'simple' phase noise ijj' of the amplifier. The 
superscripts 'prime' and 'second' are introduced in order to keep the symbols 
and ip for the overall phase fluctuation. Therefore, the phase fluctuation seen 
by the Leeson effect (Sec. [5]) is ip = tjj' + ip". After (46 1, the first phase-noise 
contribution is 

$'( s )=^l±lZlx( S ) . 
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By inspection on Fig. [24l the second phase noise contribution is 



x.it, \ s + !/ T s + V T w \ 

&'(8)=w 1 w v -—!- —X(s) . 

S + 7/T s 

Since $'(s) and i»"(s) depend deterministically on X(s), they must be added. 
Thus, the noise transfer function 



R M (:) _ *(*) («)+*"(' 



X(s) 



X(s) 



(69) 



IS 



W^Wrj s + 1/r 



1/r 



X(s). 



(70) 



The term '1' is the simple Leeson effect as introduced in Section |5j but for the 
trivial factor introduced because (70) refers to the near-dc process x f-> X 
instead of to the phase i/j if!. The term w ^ J Wv is the phase fluctuation induced 
by AM noise. 
Figures 



2G 



and 



27 



show the noise transfer function |H> C )(/)| 2 in two cases. 
The signature of the AM-PM coupling shows up in the frequency range between 
7/l and f^. In this region, the plot is parallel to that of the simple Leeson effect. 
Interestingly, the AM-PM coupling can either increase or reduce the noise in 
the region between jf L and /j> Of course, the phase-noise plot accounts for 
the slope of the near-dc process X by which the transfer function is multiplied. 
Thus for example the same signature can be seen in the l// 3 region if X is 
flicker noise. 
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Figure 27: AM-PM coupling due to parametric noise. The 'coupling' parameter 
is equal to w 7 w v /w^,. 
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Figure 28: Opto-electronic oscillator. The equivalent noise originates in both 
electronic and optical path. Gain compression takes place either in the amplifier, 
in the photodetector or in the electro-optic intensity modulator. 



9 Extended Leeson effect for the delay-line os- 
cillator 

The delay-line oscillator is a variant of the oscillator in which the resonator 
is replaced by a delay line of delay r^, so that the oscillation frequency is an 
integer multiple of To the extent of the Leeson effect, the delay line 

is equivalent to a resonator of quality factor Q — irvoTd because the slope 
d arg(/3) / 'duj is the same. Of course, longer delay gives access to lower phase 
noise and higher frequency stability, provided the delay be stable. For this 
reason the modern version of the delay-line oscillator, called OEO [T5l [TJo] and 
shown in Fig.[28j makes use of an optical fiber as the delay unit. The optical fiber 
exhibits high thermal stability (6.85x 10 _6 /K) and low loss (0.15 dB/km at 1.55 
/xm wavelength, equivalent to 0.03 dB//is), limited by the Rayleigh scattering. 
Other implementations are possible, based on a surface-wave devices and on 
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Figure 29: Noise transfer function of the opto-electronic oscillator, either for 
AM noise or PM noise. The parameter m is the mode order. Hence with 20 /is 
delay, the mode of order 2xl0 5 falls at vq — 10 GHz. 



electrical lines. 

The phase noise theory of the delay-line oscillator is widely discussed in (TJ 
Chap 5] . Here we give some key element to extend the theory to AM noise and 
to the AM-PM noise coupling. 

Since the delay line is a wide-band device, the loop can sustain oscillations 
at any frequency multiple of l/r^. A mode-selector filter is therefore necessary 
to choose one frequency by lowering the loop gain at the other frequencies. For 
this reason the feedback function of Fig. [T] is split into delay and filter, denoted 
with the subscripts d and /, respectively. It is important to understand that 
group delay of the mode selector must be orders of magnitude shorter than the 
delay of the line because the sensitivity to environment parameters is weighted 
proportionally to the delay. 

There are practical reasons to use a resonator as the mode selector. We 
assume that the delay-line attenuation is independent of frequency, moving the 
flatness defect to the resonator transfer function. Using the elementary theory of 
the Laplace transform and the material developed in Section [4j the slow- varying 
envelope representation of the feedback path is 

Mt) = 1 e-«/*/ « B/ (.) = 7 ^ (71) 

b d (t) = S(t-r d ) f+ B d (s)=e- STd (72) 
b(t) = b f (t)*b d (t) <-> B(s) =B / (s)B d (s) , (73) 

thus 
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Inserting such function in the phase-noise feedback loop we get 

H <» = rrW (75) 

and therefore 

H <»> " .-mi'-'-V, ■ <™> 

an example of which is shown in Fig. |29| 

The full extension to AM noise, as derived in Section [7] for the oscillator 
based on a simple resonator, takes cumbersome and tedious algebra. Yet at low 
frequencies, below the Leeson frequency, the asymptotic approximation of the 
delay line is a resonator of quality factor Q = nvoTd- This simplification gives 
account for the low-frequency behavior, and at least a qualitative prediction for 
the AM noise peaks. 

9.1 The impact of the laser RIN 

Let us start with the analysis of Figure [28] in open loop conditions. The path 
from the amplifier to the EOM is broken. First, we observe that the bandpass 
filter, however large, eliminates the harmonics at frequencies multiple of ujq. 
Thus the light power at the photodetector input can be described by 

P x (t) = P A [l + mcos(wot)] , (77) 

where m = J\{z) is the modulation index, J\{z) is the first-order Bessel function 
of the first kind, and z is proportional to the microwave voltage at the input of 
the intensity modulator [32 . Though the theoretical maximum is m ~ 1.164, 
in practice we get at most m = 1. The current at the photodetector input is 

J(t) = pPx(t) (78) 
= pP x [l + mcos(Lu Q t)] , (79) 

where p = is the photodetector responsivity, r\ the quantum efficiency, and 
fvjj\ the photon energy. Assuming a quantum efficiency of 0.6, the responsivity 
is of 0.75 A/W at 1.55 pro. wavelength, and of 0.64 A/W at 1.31 pm. Filtering 
out the dc component, the rms voltage across a resistor Rq at the photodetector 
output is 

v = —j= Rq p m P\ (rms voltage at the photodetector out) (80) 
v 2 

The path from the EOM to the microwave amplifier output of Fig. [28] can be 
seen as an 'amplifier,' plus a filter function. The 'amplifier' includes intensity 



modulator, photodetector, and microwave amplifier. By virtue of (80) the laser 
power affects the gain, thus the relative intensity noise (RIN) makes the gain 
fluctuate. Since in open loop 77 = Sv/v, the laser RIN induces a gain fluctuation 

v= 6 ^ s^(f) = s mN 
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Closing the loop, the results of Section [7] (AM noise) and Section [8] (AM-PM 
coupling) apply. 

Interestingly, the RIN of some lasers does not follow the polynomial law. 
Instead, slopes of —5 dB/dcc and —15 dB/dec appear in the spectrum. At 
the beginning of he process of collecting data from the literature, we suspect 
that this is typical of the distributed-feedback laser. Anyway regardless of 
the physical explanation beyond, the presence of —5 dB/dec and —15 dB/dcc 
slopes in the RIN spectrum in conjunction with the Leeson effect could explain 
the slope of —35 dB/dec and —25 dB/dec observed in the phase noise spectrum 
of some oscillators. 
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A Exotic issues 

This Appendix is not a finished work. We report on some facts intended to 
be the seed for further analysis. 

A.l AM-PM coupling in the off-resonance resonator 

It has been shown in Section [4] that the resonator operated at the exact natu- 
ral frequency responds to a phase perturbation with a decaying exponential of 
phase, with no effect on the amplitude; and that it responds to an amplitude 
perturbation with a decaying exponential of amplitude, with no effect on the 
phase. It has also been shown that cross terms appear off the resonance, for the 
resonator response is described by [b(t)] O [B(s) 



(32jH|33j). 

The AM-PM coupling inside the resonator yields naturally to the oscillator 
model depicted in Fig. 30 In this figure, the symbol s (expressed or implied) 
cannot be the Laplace complex variable because the system is nonlinear. In- 
stead, s is to be interpreted as the derivative operator 4f, which is allowed. 
Hence, the simplest approach is to derive the resonator equations using the 
Laplace formalism, and then to convert these equations into regular differential 
equations by replacing s — > jL 

The lower loop of Fig. [30| yields the Leeson effect, as described in Sec. [5] The 
upper loop models the amplitude noise as discussed in Sec. [M7j The two loops 



are coupled by the terms B aip and B va , which are nonzero when the resonator 



is off resonance. By inspection on Fig. 30 we get 



U a = £ + AB aa U a + B va ^> 



E(t) ~ S(s) 




V = 1 + Oi, 



<x v (t) ~ iys) 



X 



Uj(t) <-> Y(s) 











i 1 


^-1 — 1 — 













AM-PM 
coupling 




Figure 30: Detuning the resonator results in coupling AM to PM. 
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$ = 4- + £ $ + AB av U a 



which can be rewritten as 



1 



1 — AB n 



$ = 



1 - B„ 



1 - AB a „ 

-(J fx 



1-B„ 



Combining the above equations, we get 
1 ^ B, nry 



AB <p a B ai p 



Q l-AB aa (l-AB aa )(l-B w ) (l-AB aa ) (!-£_) Q 



$ = 



1 



1-B 



AB ai , 



-£ 



AB a <pB^p a 



<p*p 



(1- AB aa )(l-B vv ) (l-AB aa ){l-B vv ) 



hence 

U a = 



(l-AB Qa )(l-B ( 
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(81) 

The above equation is let in closed form for further analysis. The following 
formulae will be useful 



1 AB aa — 



1 - B u 



+ ±(2-A)s+(l-A)(±+n 2 ) 



s 2 + -s + 



n 2 



S 2 + is 



n 2 



n 2 



A. 2 Parametric fluctuation of the S matrix 

All over this report, the oscillator loop is analyzed as a simple block diagram in 
which the signal flows in one direction only, and there is no interaction due to 
impedances. Breaking this assumption, the amplifier and the resonator can be 
described in terms of the scatter matrix S. 

The case of the traditional microwave oscillator, where amplifier and res- 
onator are described by a 2x2 matrix, is shown in Fig. [31] The gain A and the 
resonator transfer function /3(s) are the element S21 of the respective matrix. 
Hence, the amplifier AM and PM noise as introduced in the previous Sections, 
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Figure 31: In most oscillators, the amplifier and the resonator can be described 
in terms of scatter matrix. 
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Figure 32: Scatter matrix representation of the negative-resistance (dipolar) 
oscillator. 



go in ^{^i} and ^{S^i}, respectively. The finite isolation of the amplifier is 
represented as S12 7^ 0. This effect has little importance because the isolation 
ratio of actual amplifier is high enough for the reverse signal not to circulate in 
the loop. Another effect is due to the amplifier input and output impedances, 
related to the scatter matrix by Z = y^jZq, Ti n = Su, and r out = S^- The 
amplifier input and output impedances interact with the resonator parameters. 
Thus, the fluctuations of Sn and S22 turn into frequency fluctuations. 

The quartz oscillator and other negative-resistance oscillators can also be 



described with the scatter matrix formalism (Fig. 32 1. In this case, the resonator 
degenerates into a single- value matrix. The amplifier Sn models the negative 
resistance that makes the system oscillate. Strictly, only Sn is necessary. Yet, 
in most cases the amplifier also acts as a buffer of gain S21 — ^bufj reverse gain 
(isolation) S12 = A rev , and output reflection coefficient S22 = r ou t- 

Of course, the scattering matrix formalism also apply to optics. The laser 
(oscillator) differs from the above analysis in that the laser amplifier is bidirec- 
tional and the signal may be a stationary wave. 



A.3 The Miller effect 

The Miller theorem [33] states that an impedance Z in the feedback path of an 
amplifier of gain A can be replaced by two impedances, 



I- A 



(input) 



and 



1 - l/A 



(output) 
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Figure 33: Miller effect. 



connected at the amplifier input output, respectively (Fig. 33 1. For our pur- 
poses, the left-hand side of Fig. [33] is formally equivalent to the oscillator loop, 
for we can identify Z as the resonator in the feedback path, and A as the sus- 
taining amplifier. 

Unfortunately the Miller theorem cannot be inverted in the general case 
because it would be necessary to collapse three degrees of freedom (A, Z\ and 
Z2) into two degrees of freedom (A and Z). The parameters of the specific 
circuit are needed to get Z from Z\ and Z2. Nonetheless, the Miller theorem 
provides evidence that the gain fluctuations affect the impedances of the whole 
circuit, and that a fluctuating impedance at the amplifier input or output can 
be turned into a fluctuating impedance in parallel to Z, hence to the resonator. 
In turn, the oscillator frequency fluctuates. 
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